Paired metagenomic and chemical evaluation of aflatoxin-contaminated dog kibble.

Introduction: Identification of chemical toxins from complex or highly processed foods can present 'needle in the haystack' challenges for chemists. Metagenomic data can be used to guide chemical toxicity evaluations by providing DNA-based description of the wholistic composition (eukaryotic, bacterial, protozoal, viral, and antimicrobial resistance) of foods suspected to harbor toxins, allergens, or pathogens. This type of information can focus chemistry-based diagnostics, improve hazard characterization and risk assessment, and address data gaps. Additionally, there is increasing recognition that simultaneously co-occurring mycotoxins, either from single or multiple species, can impact dietary toxicity exposure. Metagenomic data provides a way to address data gaps related to co-occurrence of multiple fungal species. Methods: Paired metagenomic and chemical data were used to evaluate aflatoxin-contaminated kibble with known levels of specific mycotoxins. Kibble was ground to a fine powder for both chemical and molecular analyses. Chemical analyses were performed with Liquid Chromatography Mass Spectrometry (LCMS) and according to the AOAC Official method 2005.08: Aflatoxins in Corn, Raw Peanuts, and Peanut Butter using Liquid Chromatography with Post-Column Photochemical Derivatization. Metagenomes were created from DNA extracted from ground kibble and sequenced on an Illumina NextSeq 2000 with an average sequence depth of 180 million reads per replicate. Results and discussion: Metagenomic data demonstrated that the abundance of DNA from putative aflatoxigenic Aspergillus spp. correlated with the levels of aflatoxin quantified by LCMS. Metagenomic data also identified an expansive range of co-occurring fungal taxa which may produce additional mycotoxins. DNA data paired with chemical data provides a novel modality to address current data gaps surrounding dietary mycotoxin exposure, toxigenic fungal taxonomy, and mycotoxins of emerging concern.

The spread of pESI-mediated extended-spectrum cephalosporin resistance in Salmonella serovars-Infantis, Senftenberg, and Alachua isolated from food animal sources in the United States.

The goal of this study is to investigate the origin, prevalence, and evolution of the pESI megaplasmid in Salmonella isolated from animals, foods, and humans. We queried 510,097 Salmonella genomes under the National Center for Biotechnology Information (NCBI) Pathogen Detection (PD) database for the presence of potential sequences containing the pESI plasmid in animal, food, and environmental sources. The presence of the pESI megaplasmid was confirmed by using seven plasmid-specific markers (rdA, pilL, SogS, TrbA, ipf, ipr2 and IncFIB(pN55391)). The plasmid and chromosome phylogeny of these isolates was inferred from single nucleotide polymorphisms (SNPs). Our search resolved six Salmonella clusters carrying the pESI plasmid. Four were emergent Salmonella Infantis clusters, and one each belonged to serovar Senftenberg and Alachua. The Infantis cluster with a pESI plasmid carrying blaCTX-M-65 gene was the biggest of the four emergent Infantis clusters, with over 10,000 isolates. This cluster was first detected in South America and has since spread widely in United States. Over time the composition of pESI in United States has changed with the average number of resistance genes showing a decrease from 9 in 2014 to 5 in 2022, resulting from changes in gene content in two integrons present in the plasmid. A recent and emerging cluster of Senftenberg, which carries the blaCTX-M-65 gene and is primarily associated with turkey sources, was the second largest in the United States. SNP analysis showed that this cluster likely originated in North Carolina with the recent acquisition of the pESI plasmid. A single Alachua isolate from turkey was also found to carry the pESI plasmid containing blaCTX-M-65 gene. The study of the pESI plasmid, its evolution and mechanism of spread can help us in developing appropriate strategies for the prevention and further spread of this multi-drug resistant plasmid in Salmonella in poultry and humans.

Metagenomic survey of antimicrobial resistance (AMR) in Maryland surface waters differentiated by high and low human impact.

Here, we examine surface waters as a modality to better understand baseline antimicrobial resistance (AMR) across the environment to supplement existing AMR monitoring in pathogens associated with humans, foods, and animals. Data from metagenomic and quasimetagenomic (shotgun sequenced enrichments) are used to describe AMR in Maryland surface waters from high and low human impact classifications.

Fecal microbiomes of laboratory beagles receiving antiparasitic formulations in an experimental setting.

Here, we describe the fecal microbiome of laboratory beagles in a non-invasive experiment designed to contrast in vivo versus in vitro bioequivalence in response to antiparasitic drug administration. The experiment provided a unique opportunity to evaluate metagenomic profiles of canine feces before and after anti-parasitic drug exposure.

Genomic analysis of azithromycin-resistant Salmonella from food animals at slaughter and processing, and retail meats, 2011-2021, United States.

IMPORTANCE: Macrolides of different ring sizes are critically important antimicrobials for human medicine and veterinary medicine, though the widely used 15-membered ring azithromycin in humans is not approved for use in veterinary medicine. We document here the emergence of azithromycin-resistant Salmonella among the NARMS culture collections between 2011 and 2021 in food animals and retail meats, some with co-resistance to ceftriaxone or decreased susceptibility to ciprofloxacin. We also provide insights into the underlying genetic mechanisms and genomic contexts, including the first report of a novel combination of azithromycin resistance determinants and the characterization of multidrug-resistant plasmids. Further, we highlight the emergence of a multidrug-resistant Salmonella Newport clone in food animals (mainly cattle) with both azithromycin resistance and decreased susceptibility to ciprofloxacin. These findings contribute to a better understating of azithromycin resistance mechanisms in Salmonella and warrant further investigations on the drivers behind the emergence of resistant clones.

Application of quasimetagenomics methods to define microbial diversity and subtype Listeria monocytogenes in dairy and seafood production facilities.

IMPORTANCE: In developed countries, the human diet is predominated by food commodities, which have been manufactured, processed, and stored in a food production facility. Little is known about the application of metagenomic sequencing approaches for detecting foodborne pathogens, such as L. monocytogenes, and characterizing microbial diversity in food production ecosystems. In this work, we investigated the utility of 16S rRNA amplicon and quasimetagenomic sequencing for the taxonomic and phylogenetic classification of Listeria culture enrichments of environmental swabs collected from dairy and seafood production facilities. We demonstrated that single-nucleotide polymorphism (SNP) analyses of L. monocytogenes metagenome-assembled genomes (MAGs) from quasimetagenomic data sets can achieve similar resolution as culture isolate whole-genome sequencing. To further understand the impact of genome coverage on MAG SNP cluster resolution, an in silico downsampling approach was employed to reduce the percentage of target pathogen sequence reads, providing an initial estimate of required MAG coverage for subtyping resolution of L. monocytogenes.

Genotypic and Phenotypic Characterization of Antimicrobial and Heavy Metal Tolerance in Salmonella enterica and Escherichia coli Isolates from Swine Feed Mills.

Antimicrobials and heavy metals are commonly used in the animal feed industry. The role of in-feed antimicrobials on the evolution and persistence of resistance in enteric bacteria is not well described. Whole-Genome Sequencing (WGS) is widely used for genetic characterizations of bacterial isolates, including antimicrobial resistance, heavy metal tolerance, virulence factors, and relatedness to other sequenced isolates. The goals of this study were to i) use WGS to characterize Salmonella enterica (n = 33) and Escherichia coli (n = 30) isolated from swine feed and feed mill environments; and ii) investigate their genotypic and phenotypic antimicrobial and heavy metal tolerance. Salmonella isolates belonged to 10 serovars, the most common being Cubana, Senftenberg, and Tennessee. E. coli isolates were grouped into 22 O groups. Phenotypic resistance to at least one antimicrobial was observed in 19 Salmonella (57.6%) and 17 E. coli (56.7%) isolates, whereas multidrug resistance (resistant to ≥3 antimicrobial classes) was observed in four Salmonella (12%) and two E. coli (7%) isolates. Antimicrobial resistance genes were identified in 17 Salmonella (51%) and 29 E. coli (97%), with 11 and 29 isolates possessing genes conferring resistance to multiple antimicrobial classes. Phenotypically, 53% Salmonella and 58% E. coli presented resistance to copper and arsenic. All isolates that possessed the copper resistance operon were resistant to the highest concentration tested (40 mM). Heavy metal tolerance genes to copper and silver were present in 26 Salmonella isolates. Our study showed a strong agreement between predicted and measured resistances when comparing genotypic and phenotypic data for antimicrobial resistance, with an overall concordance of 99% and 98.3% for Salmonella and E. coli, respectively.

Lociq provides a loci-seeking approach for enhanced plasmid subtyping and structural characterization.

Antimicrobial resistance (AMR) monitoring for public health is relying more on whole genome sequencing to characterize and compare resistant strains. This requires new approaches to describe and track AMR that take full advantage of the detailed data provided by genomic technologies. The plasmid-mediated transfer of AMR genes is a primary concern for AMR monitoring because plasmid rearrangement events can integrate new AMR genes into the plasmid backbone or promote hybridization of multiple plasmids. To better monitor plasmid evolution and dissemination, we developed the Lociq subtyping method to classify plasmids by variations in the sequence and arrangement of core plasmid genetic elements. Subtyping with Lociq provides an alpha-numeric nomenclature that can be used to denominate plasmid population diversity and characterize the relevant features of individual plasmids. Here we demonstrate how Lociq generates typing schema to track and characterize the origin, evolution and epidemiology of multidrug resistant plasmids.

Geography shapes the genomics and antimicrobial resistance of Salmonella enterica Serovar Enteritidis isolated from humans.

Multidrug-resistant (MDR) Salmonella has been a long-standing challenge in public health and food safety. The prevalence of MDR S. Enteritidis, especially isolated from humans, in China is significantly higher than those from the U.S. and other countries. A dataset of 197 S. Enteritidis genomes, including 16 sequenced clinical isolates from China and 181 downloaded genomes of human isolates from the U.S., Europe, and Africa, was analyzed for genomic diversity, virulence potential, and antimicrobial resistance (AMR). Phylogenomic analyses identified four major well-supported clades (I-IV). While AMR genotype in the majority of isolates in clades I and IV displayed as pan-susceptible, 81.8% (9/11) and 22.4% (13/58) of isolates in clades III and II were MDR, respectively. It is noted that 77% (10/13) of MDR isolates in clade II were from China. The most common antimicrobial resistance genes (ARGs) carried by the Chinese isolates were aph(3')-IIa, blaCTX-M-55, and blaTEM-1B, whereas blaTEM-1B, sul1, sul2, drfA7, aph(3")-Ib/strA, and aph(6)-Id/strB were most often identified in those from Africa (clade III). Among the 14 plasmid types identified, IncX1 and IncFII(pHN7A8) were found exclusively in the Chinese MDR isolates, while IncQ1 was highly associated with the African MDR isolates. The spvRABCD virulence operon was present in 94.9% (187/197) of isolates tested and was highly associated with both the IncF (IncFII and IncFIB) plasmids. In addition, phylogenetic differences in distribution of Salmonella pathogenicity islands (SPIs), prophages and other accessory genes were also noted. Taken together, these findings provide new insights into the molecular mechanisms underpinning diversification of MDR S. Enteritidis.

PRAWNS: compact pan-genomic features for whole-genome population genomics.

MOTIVATION: Scientists seeking to understand the genomic basis of bacterial phenotypes, such as antibiotic resistance, today have access to an unprecedented number of complete and nearly complete genomes. Making sense of these data requires computational tools able to perform multiple-genome comparisons efficiently, yet currently available tools cannot scale beyond several tens of genomes. RESULTS: We describe PRAWNS, an efficient and scalable tool for multiple-genome analysis. PRAWNS defines a concise set of genomic features (metablocks), as well as pairwise relationships between them, which can be used as a basis for large-scale genotype-phenotype association studies. We demonstrate the effectiveness of PRAWNS by identifying genomic regions associated with antibiotic resistance in Acinetobacter baumannii. AVAILABILITY AND IMPLEMENTATION: PRAWNS is implemented in C++ and Python3, licensed under the GPLv3 license, and freely downloadable from GitHub (https://github.com/KiranJavkar/PRAWNS.git). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Listeria monocytogenes Illness and Deaths Associated With Ongoing Contamination of a Multiregional Brand of Ice Cream Products, United States, 2010-2015.

BACKGROUND: Frozen foods have rarely been linked to Listeria monocytogenes illness. We describe an outbreak investigation prompted by both hospital clustering of illnesses and product testing. METHODS: We identified outbreak-associated listeriosis cases using whole-genome sequencing (WGS), product testing results, and epidemiologic linkage to cases in the same Kansas hospital. We reviewed hospital medical and dietary records, product invoices, and molecular subtyping results. Federal and state officials tested product and environmental samples for L. monocytogenes. RESULTS: Kansas officials were investigating 5 cases of listeriosis at a single hospital when, simultaneously, unrelated sampling for a study in South Carolina identified L. monocytogenes in Company A ice cream products made in Texas. Isolates from 4 patients and Company A products were closely related by WGS, and the 4 patients with known exposures had consumed milkshakes made with Company A ice cream while hospitalized. Further testing identified L. monocytogenes in ice cream produced in a second Company A production facility in Oklahoma; these isolates were closely related by WGS to those from 5 patients in 3 other states. These 10 illnesses, involving 3 deaths, occurred from 2010 through 2015. Company A ultimately recalled all products. CONCLUSIONS: In this US outbreak of listeriosis linked to a widely distributed brand of ice cream, WGS and product sampling helped link cases spanning 5 years to 2 production facilities, indicating longstanding contamination. Comprehensive sanitation controls and environmental and product testing for L. monocytogenes with regulatory oversight should be implemented for ice cream production.

Toward the Adoption of Loop-Mediated Isothermal Amplification for Salmonella Screening at the National Antimicrobial Resistance Monitoring System's Retail Meat Sites.

The National Antimicrobial Resistance Monitoring System (NARMS) is a One Health program in the United States that collects data on antimicrobial resistance in enteric bacteria from humans, animals, and the environment. Salmonella is a major pathogen tracked by the NARMS retail meat arm but currently lacks a uniform screening method. We evaluated a loop-mediated isothermal amplification (LAMP) assay for the rapid screening of Salmonella from 69 NARMS retail meat and poultry samples. All samples were processed side by side for culture isolation using two protocols, one from NARMS and the other one described in the U.S. Food and Drug Administration's Bacteriological Analytical Manual (BAM). Overall, 10 (14.5%) samples screened positive by the Salmonella LAMP assay. Of those, six were culture-confirmed by the NARMS protocol and six by the BAM method with overlap on four samples. No Salmonella isolates were recovered from samples that screened negative with LAMP. These results suggested 100% sensitivity for LAMP in reference to culture. Antimicrobial susceptibility testing and whole-genome sequencing analysis confirmed identities of these isolates. Using the BAM protocol, all Salmonella isolates were recovered from samples undergoing Rappaport-Vassiliadis medium selective enrichment and presumptive colonies (n = 130) were dominated by Hafnia alvei (44.6%), Proteus mirabilis (22.3%), and Morganella morganii (9.9%) based on matrix-assisted laser desorption/ionization time-of-flight mass spectrometry. This method comparison study clearly demonstrated the benefit of a rapid, robust, and highly sensitive molecular screening method in streamlining the laboratory workflow. Fourteen NARMS retail meat sites further verified the performance of this assay using a portion of their routine samples, reporting an overall specificity of 98.8% and sensitivity of 90%. As of July 2022, the vast majority of NARMS retail meat sites have adopted the Salmonella LAMP assay for rapid screening of Salmonella in all samples.

Prevalence of Antimicrobial Resistance in Select Bacteria From Retail Seafood-United States, 2019.

In 2019, the United States National Antimicrobial Resistance Monitoring System (NARMS) surveyed raw salmon, shrimp, and tilapia from retail grocery outlets in eight states to assess the prevalence of bacterial contamination and antimicrobial resistance (AMR) in the isolates. Prevalence of the targeted bacterial genera ranged among the commodities: Salmonella (0%-0.4%), Aeromonas (19%-26%), Vibrio (7%-43%), Pseudomonas aeruginosa (0.8%-2.3%), Staphylococcus (23%-30%), and Enterococcus (39%-66%). Shrimp had the highest odds (OR: 2.8, CI: 2.0-3.9) of being contaminated with at least one species of these bacteria, as were seafood sourced from Asia vs. North America (OR: 2.7; CI: 1.8-4.7) and Latin America and the Caribbean vs. North America (OR: 1.6; CI: 1.1-2.3) and seafood sold at the counter vs. sold frozen (OR: 2.1; CI: 1.6-2.9). Isolates exhibited pan-susceptibility (Salmonella and P. aeruginosa) or low prevalence of resistance (<10%) to most antimicrobials tested, with few exceptions. Seafood marketed as farm-raised had lower odds of contamination with antimicrobial resistant bacteria compared to wild-caught seafood (OR: 0.4, CI: 0.2-0.7). Antimicrobial resistance genes (ARGs) were detected for various classes of medically important antimicrobials. Clinically relevant ARGs included carbapenemases (bla IMI-2, bla NDM-1) and extended spectrum ß-lactamases (ESBLs; bla CTX-M-55). This population-scale study of AMR in seafood sold in the United States provided the basis for NARMS seafood monitoring, which began in 2020.

Use of Large-Scale Genomics to Identify the Role of Animals and Foods as Potential Sources of Extraintestinal Pathogenic Escherichia coli That Cause Human Illness

Extraintestinal pathogenic Escherichia coli (ExPEC) cause urinary tract and potentially life-threatening invasive infections. Unfortunately, the origins of ExPEC are not always clear. We used genomic data of E. coli isolates from five U.S. government organizations to evaluate potential sources of ExPEC infections. Virulence gene analysis of 38,032 isolates from human, food animal, retail meat, and companion animals classified the subset of 8142 non-diarrheagenic isolates into 40 virulence groups. Groups were identified as low, medium, and high relative risk of containing ExPEC strains, based on the proportion of isolates recovered from humans. Medium and high relative risk groups showed a greater representation of sequence types associated with human disease, including ST-131. Over 90% of food source isolates belonged to low relative risk groups, while >60% of companion animal isolates belonged to medium or high relative risk groups. Additionally, 18 of the 26 most prevalent antimicrobial resistance determinants were more common in high relative risk groups. The associations between antimicrobial resistance and virulence potentially limit treatment options for human ExPEC infections. This study demonstrates the power of large-scale genomics to assess potential sources of ExPEC strains and highlights the importance of a One Health approach to identify and manage these human pathogens.

Multi-laboratory evaluation of the Illumina iSeq platform for whole genome sequencing of Salmonella, Escherichia coli and Listeria.

There is a growing need for public health and veterinary laboratories to perform whole genome sequencing (WGS) for monitoring antimicrobial resistance (AMR) and protecting the safety of people and animals. With the availability of smaller and more affordable sequencing platforms coupled with well-defined bioinformatic protocols, the technological capability to incorporate this technique for real-time surveillance and genomic epidemiology has greatly expanded. There is a need, however, to ensure that data are of high quality. The goal of this study was to assess the utility of a small benchtop sequencing platform using a multi-laboratory verification approach. Thirteen laboratories were provided the same equipment, reagents, protocols and bacterial reference strains. The Illumina DNA Prep and Nextera XT library preparation kits were compared, and 2×150 bp iSeq i100 chemistry was used for sequencing. Analyses comparing the sequences produced from this study with closed genomes from the provided strains were performed using open-source programs. A detailed, step-by-step protocol is publicly available via protocols.io (https://www.protocols.io/view/iseq-bacterial-wgs-protocol-bij8kcrw). The throughput for this method is approximately 4-6 bacterial isolates per sequencing run (20-26 Mb total load). The Illumina DNA Prep library preparation kit produced high-quality assemblies and nearly complete AMR gene annotations. The Prep method produced more consistent coverage compared to XT, and when coverage benchmarks were met, nearly all AMR, virulence and subtyping gene targets were correctly identified. Because it reduces the technical and financial barriers to generating WGS data, the iSeq platform is a viable option for small laboratories interested in genomic surveillance of microbial pathogens.

Long-Read Sequencing Reveals Evolution and Acquisition of Antimicrobial Resistance and Virulence Genes in Salmonella enterica.

Salmonella enterica is a significant and phylogenetically diverse zoonotic pathogen. To understand its genomic heterogeneity and antimicrobial resistance, we performed long-read sequencing on Salmonella isolated from retail meats and food animals. A collection of 134 multidrug-resistant isolates belonging to 33 serotypes were subjected to PacBio sequencing. One major locus of diversity among these isolates was the presence and orientation of Salmonella pathogenic islands (SPI), which varied across different serotypes but were largely conserved within individual serotypes. We also identified insertion of an IncQ resistance plasmid into the chromosome of fourteen strains of serotype I 4,[5],12:i:- and the Salmonella genomic island 1 (SGI-1) in five serotypes. The presence of various SPIs, SGI-1 and integrated plasmids contributed significantly to the genomic variability and resulted in chromosomal resistance in 55.2% (74/134) of the study isolates. A total of 93.3% (125/134) of isolates carried at least one plasmid, with isolates carrying up to seven plasmids. We closed 233 plasmid sequences of thirteen replicon types, along with twelve hybrid plasmids. Some associations between Salmonella isolate source, serotype, and plasmid type were seen. For instance, IncX plasmids were more common in serotype Kentucky from retail chicken. Plasmids IncC and IncHI had on average more than five antimicrobial resistance genes, whereas in IncX, it was less than one per plasmid. Overall, 60% of multidrug resistance (MDR) strains that carried >3 AMR genes also carried >3 heavy metal resistance genes, raising the possibility of co-selection of antimicrobial resistance in the presence of heavy metals. We also found nine isolates representing four serotypes that carried virulence plasmids with the spv operon. Together, these data demonstrate the power of long-read sequencing to reveal genomic arrangements and integrated plasmids with a high level of resolution for tracking and comparing resistant strains from different sources. Additionally, the findings from this study will help expand the reference set of closed Salmonella genomes that can be used to improve genome assembly from short-read data commonly used in One Health antimicrobial resistance surveillance.

Whole-Genome Assessment of Clinical Acinetobacter baumannii Isolates Uncovers Potentially Novel Factors Influencing Carbapenem Resistance.

Carbapenems-one of the important last-line antibiotics for the treatment of gram-negative infections-are becoming ineffective for treating Acinetobacter baumannii infections. Studies have identified multiple genes (and mechanisms) responsible for carbapenem resistance. In some A. baumannii strains, the presence/absence of putative resistance genes is not consistent with their resistance phenotype-indicating the genomic factors underlying carbapenem resistance in A. baumannii are not fully understood. Here, we describe a large-scale whole-genome genotype-phenotype association study with 349 A. baumannii isolates that extends beyond the presence/absence of individual antimicrobial resistance genes and includes the genomic positions and pairwise interactions of genes. Ten known resistance genes exhibited statistically significant associations with resistance to imipenem, a type of carbapenem: blaOXA-23, qacEdelta1, sul1, mphE, msrE, ant(3")-II, aacC1, yafP, aphA6, and xerD. A review of the strains without any of these 10 genes uncovered a clade of isolates with diverse imipenem resistance phenotypes. Finer resolution evaluation of this clade revealed the presence of a 38.6 kbp conserved chromosomal region found exclusively in imipenem-susceptible isolates. This region appears to host several HTH-type DNA binding transcriptional regulators and transporter genes. Imipenem-susceptible isolates from this clade also carried two mutually exclusive plasmids that contain genes previously known to be specific to imipenem-susceptible isolates. Our analysis demonstrates the utility of using whole genomes for genotype-phenotype correlations in the context of antibiotic resistance and provides several new hypotheses for future research.

Core Genome MLST for Source Attribution of Campylobacter coli.

Campylobacter species are among the leading foodborne bacterial agents of human diarrheal illness. The majority of campylobacteriosis has been attributed to Campylobacter jejuni (85% or more), followed by Campylobacter coli (5-10%). The distribution of C. jejuni and C. coli varies by host organism, indicating that the contribution to human infection may differ between isolation sources. To address the relative contribution of each source to C. coli infections in humans, core genome multilocus sequence type with a 200-allele difference scheme (cgMLST200) was used to determine cgMLST type for 3,432 C. coli isolated from food animals (n = 2,613), retail poultry meats (n = 389), human clinical settings (n = 285), and environmental sources (n = 145). Source attribution was determined by analyzing the core genome with a minimal multilocus distance methodology (MMD). Using MMD, a higher proportion of the clinical C. coli population was attributed to poultry (49.6%) and environmental (20.9%) sources than from cattle (9.8%) and swine (3.2%). Within the population of C. coli clinical isolates, 70% of the isolates that were attributed to non-cecal retail poultry, dairy cattle, beef cattle and environmental waters came from two cgMLST200 groups from each source. The most common antibiotic resistance genes among all C. coli were tetO (65.6%), bla OXA - 193 (54.2%), aph(3')-IIIa (23.5%), and aadE-Cc (20.1%). Of the antibiotic resistance determinants, only one gene was isolated from a single source: bla OXA - 61 was only isolated from retail poultry. Within cgMLST200 groups, 17/17 cgMLST200-435 and 89/92 cgMLST200-707 isolates encoded for aph(3')-VIIa and 16/16 cgMLST200-319 harbored aph(2')-If genes. Distribution of bla OXA alleles showed 49/50 cgMLST200-5 isolates contained bla OXA - 498 while bla OXA - 460 was present in 37/38 cgMLST200-650 isolates. The cgMLST200-514 group revealed both ant(6)-Ia and sat4 resistance genes in 23/23 and 22/23 isolates, respectively. Also, cgMLST200-266 and cgMLST200-84 had GyrAT86I mutation with 16/16 (100%) and 14/15 (93.3%), respectively. These findings illustrate how cgMLST and MMD methods can be used to evaluate the relative contribution of known sources of C. coli to the human burden of campylobacteriosis and how cgMLST typing can be used as an indicator of antimicrobial resistance in C. coli.

Evaluating the accuracy of Listeria monocytogenes assemblies from quasimetagenomic samples using long and short reads.

BACKGROUND: Whole genome sequencing of cultured pathogens is the state of the art public health response for the bioinformatic source tracking of illness outbreaks. Quasimetagenomics can substantially reduce the amount of culturing needed before a high quality genome can be recovered. Highly accurate short read data is analyzed for single nucleotide polymorphisms and multi-locus sequence types to differentiate strains but cannot span many genomic repeats, resulting in highly fragmented assemblies. Long reads can span repeats, resulting in much more contiguous assemblies, but have lower accuracy than short reads. RESULTS: We evaluated the accuracy of Listeria monocytogenes assemblies from enrichments (quasimetagenomes) of naturally-contaminated ice cream using long read (Oxford Nanopore) and short read (Illumina) sequencing data. Accuracy of ten assembly approaches, over a range of sequencing depths, was evaluated by comparing sequence similarity of genes in assemblies to a complete reference genome. Long read assemblies reconstructed a circularized genome as well as a 71 kbp plasmid after 24 h of enrichment; however, high error rates prevented high fidelity gene assembly, even at 150X depth of coverage. Short read assemblies accurately reconstructed the core genes after 28 h of enrichment but produced highly fragmented genomes. Hybrid approaches demonstrated promising results but had biases based upon the initial assembly strategy. Short read assemblies scaffolded with long reads accurately assembled the core genes after just 24 h of enrichment, but were highly fragmented. Long read assemblies polished with short reads reconstructed a circularized genome and plasmid and assembled all the genes after 24 h enrichment but with less fidelity for the core genes than the short read assemblies. CONCLUSION: The integration of long and short read sequencing of quasimetagenomes expedited the reconstruction of a high quality pathogen genome compared to either platform alone. A new and more complete level of information about genome structure, gene order and mobile elements can be added to the public health response by incorporating long read analyses with the standard short read WGS outbreak response.